The besy1() function returns the y1st Bessel function of its argument, i.e it returns the first order Bessel function of the second kind. besy1() expects its argument to be in radians.